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I.  INTRODUCTION 


Sonic  and  ultrasonic  techniques  are  available  for  determining  the 
acoustic  properties  of  materials  b  '  making  direct  laboratory  measure¬ 
ments  of  the  phase  velocity  and  attenuation  of  harmonic  waves.  Numerous 
authors  have  reported  such  measurements  on  saturated  sediments.  Applied 
Research  Laboratories,  The  University  of  Texas  at  Austin  (ARL:UT),  has 
carried  out  many  measurements  of  this  type  and  has  pioneered  the 
development  of  transducers  for  the  generation  and  sensing  of  shear 
waves. 

There  is  currently  controversy  over  the  frequency  dependence  of  the 

shear  and  compressional  attenuations  (a°d  hence,  on  the  basis  of 

causality  arguments,  their  velocities)  of  saturated  marine  sediments. 

This  frequency  dependence  is  an  important  basic  property  of  sand 

sediments  that  could  be  related  to  the  unexplained  frequency  dependence 

of  low  frequency  propagation  in  shallow  water. ®  Hamilton,^  using  a 

viscoelastic  framework,  carried  out  an  extensive  examination  of  reported 

attenuation  values  for  compressional  waves  in  marine  sediments  and 

concluded  that  these  values  were  consistent  with  a  linear  dependence  on 

frequency.  Extensive  analytical  predictions  of  the  phase  velocity  and 

attenuation  of  acoustic  waves  in  saturated  sediments  have  been  presented 

by  Stoll  and  Bryan8  and  by  Stoll, who  based  their  work  on  the 
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classical  mixture  theory  of  Biot.  ’  The  theory  of  Stoll  and  Bryan  was 

4 

further  refined  by  Hovem  and  Ingram,  and  by  Bedford,  Costley,  and 

14 

Stern.  These  theories,  based  on  the  microstructure  of  the  sediment, 
predict  a  more  complicated  frequency  dependence  than  that  suggested  by 
Hamilton.  In  the  absence  of  direct  measurements  at  low  frequencies, 
experiments  that  test  the  validity  of  the  Biot  approach  are  needed  to 
determine  the  correct  method  to  use  in  extrapolating  measurements  at 
higher  frequencies  to  lower  frequencies. 

15 

ARL:UT  recently  introduced  an  important  innovation  in  the  study 
of  acoustic  properties  of  saturated  sediments.  A  mixture  of  water  and 
glycerine  was  used  as  the  saturating  fluid;  by  varying  the  ratio  of 


glycerine  to  water,  the  viscosity,  density,  and  bulk  modulus  of  the 
saturant  could  be  varied  in  a  controlled  manner.  Since  these  parameters 
appear  explicitly  in  the  Biot-Stoll  equations,  this  procedure  permits  a 
definitive  test  of  the  theory.  Furthermore,  these  measurements  can  be 
used  in  a  more  complete  theoretical  analysis  of  the  coefficients  in  the 
equations.  Preliminary  measurements^  indicated  that  the  Biot-Stoll 
predictions  for  the  dependence  on  viscosity  may  not  be  correct. 
However,  due  to  large  experimental  uncertainty,  no  definitive 
conclusions  could  be  drawn. 

This  report  summarizes  further  laboratory  measurements  at  ARL : UT 
during  the  period  1  September  1983  -  31  August  1984.  The  focus  of  this 
research  was  to  improve  the  accuracy  of  velocity  and  attenuation 
measurements  so  that  a  definitive  test  of  the  Biot-Stoll  theory  could  be 
made.  The  major  accomplishments  during  this  period  were  refining  the 
design  of  transducers,  improving  the  procedures  for  making  accurate 
laboratory  measurements  of  the  velocity  and  attenuation  of  shear  and 
compressional  waves  in  saturated  sands,  and  locating  the  sources  of 
experimental  error.  Preliminary  measurements  of  compressional  wave 
parameters  were  also  made  as  the  first  step  in  testing  the  Biot-Stoll 
theory.  This  work  was  carried  out  by  D.  Costley  under  the  direction  of 
Dr.  A.  Bedford,  Department  of  Aerospace  Engineering  and  Engineering 
Mechanics,  The  University  of  Texas  at  Austin,  and  Dr.  Paul  Vidmar, 
Applied  Research  Laboratories,  The  University  of  Texas  at  Austin.  This 
work  has  been  submitted  by  Mr.  Costley  as  part  of  a  Master's  thesis  in 
Aerospace  and  Engineering  Mechanics  at  The  University  of  Texas  at 
Austin. 

The  remainder  of  this  report  is  organized  as  follows.  Section  II 
reviews  the  measurement  procedure  and  preliminary  results  (forming  part 
of  Mr.  Costley's  Master's  thesis).  Section  III  discusses  the  progress 
made  and  suggests  directions  for  future  research. 


II.  EXPERIMENTAL  PROCEDURE  AND  PRELIMINARY  RESULTS 


A.  General  Procedure 


The  objective  of  this  experiment  is  to  determine  the  phase  velocity 
and  attenuation  of  harmonic  compressional  waves  in  a  fluid  saturated 
sediment  as  a  function  of  viscosity  of  the  pore  fluid.  These  results 
will  be  compared  with  those  predicted  by  the  Biot-Stoll  theory.  Since 
the  parameters  viscosity,  bulk  modulus,  and  density  appear  in  the  Biot- 
Stoll  expression  for  the  wave  number,  these  measur .nents  are  expected  to 
be  a  critical  test  of  the  theory. 

The  main  focus  of  the  work  described  here  has  been  to  improve  the 
techniques  developed  by  previous  investigators^’^  for  measuring  the 
attenuation  and  phase  velocity  of  compressional  waves  and  to  locate  the 
causes  of  experimental  error.  The  most  significant  development  over 
earlier  work  was  the  redesign  of  the  compressional  wave  transducers. 
With  this  improvement,  spherical  spreading  was  established  and  measured, 
thus  reducing  the  error  caused  by  the  directivity  of  the  transducers. 
Two  identical  transducers  were  built,  one  serving  as  a  transmitter  and 
the  other  as  a  receiver.  The  transmitter  converts  an  electrical  signal, 
which  is  generated  in  an  oscillator,  into  an  acoustic  signal,  which 
propagates  through  the  sediment.  The  receiver  detects  this  acoustic 
signal  and  reconverts  it  into  an  electrical  signal  which  is  displayed  on 
an  osci 1 loscope. 

The  central  idea  in  our  measurement  procedure  is  to  compare  the 
time  delay  and  amplitude  of  the  output  voltage  signal  at  one  transducer 
separation  with  the  same  measurements  at  other  separations.  From  the 
time  delay  measurements  at  different  separations,  the  phase  velocity  is 
determined.  Similarly,  by  comparing  the  different  amplitude  measure¬ 
ments,  the  attenuation  can  be  calculated.  It  is  assumed  that  the 
coupling  between  the  transducer  and  the  sediment  is  the  same  at  each 
separation.  This  assumption  allows  us  to  relate  the  decay  in  amplitude 
only  to  the  increase  in  transducer  separation. 


Initially,  the  transducers  are  placed  approximately  20  cm  apart. 
The  time  delay  from  the  beginning  of  the  triggered  pulse  to  the  fourth 
positive  peak  of  the  received  signal  is  then  recorded  along  with  the 
peak-to-peak  amplitude  of  the  fourth  cycle  of  the  received  signal.  The 
fourth  cycle  is  used  so  that  the  measurements  are  made  as  close  to  the 
steady  state  part  of  the  signal  as  possible.  The  signal  is  then 
monitored  for  at  least  four  hours  to  ensure  that  there  are  no  variations 
in  signal  amplitude  due  to  changes  in  the  coupling  between  the  sediment 
and  the  transducers  or  due  to  any  settling  of  the  sediment.  Four  to  six 
readings  are  taken  over  this  period,  and  the  slight  variations  about  the 
stabilized  value  are  averaged.  The  receiver  is  then  moved  2.0  cm  closer 
to  the  transmitter  and  the  amplitude  and  time  delay  measurements  are 
repeated.  These  measurements  are  repeated  at  2.0  cm  intervals  until  the 
separation  between  the  transducers  is  approximately  6.0  cm.  Using  this 
set  of  measurements,  the  phase  velocity  and  attenuation  are  determined 
by  finding  the  least  squares  fit  to  the  data. 

The  viscosity  of  the  fluid  is  then  increased  by  adding  glycerine  to 
the  mixture.  The  entire  set  of  measurements  is  repeated  and  the  phase 
velocity  and  attenuation  are  determined  for  a  new  value  of  viscosity. 
This  procedure  is  repeated  several  times,  each  time  with  a  new  value  of 
viscosity.  In  this  way,  plots  of  phase  velocity  and  attenuation  versus 
viscosity  can  be  made  and  then  compared  with  predictions  of  the  theory. 

3.  Instrumentation 


A  schematic  of  the  laboratory  equipment  is  shown  in  Fig.  1.  A 
continuous  tone  of  105.0  kHz  and  10.0  V  peak-to-peak  amplitude  is 

generated  in  the  oscillator.  The  output  is  sent  to  a  frequency  counter, 
a  tone  burst  generator,  and  an  oscilloscope.  The  tone  burst  generator 
passes  a  four-cycle  pulse  of  105.0  kHz  to  the  transmitting  transducer. 
The  signal  is  then  cut  off  for  several  milliseconds  so  that  the 

reflections  of  the  signal  off  the  sides  of  the  container  and  the 

sediment-water  interface  will  die  down  before  a  new  four-cycle  pulse  is 
sent  to  the  transmitter.  Any  reflections  not  suppressed  will  interfere 
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with  the  signal  that  follows  the  direct  path  between  the  transducers, 
and  thus  produce  an  error  in  the  signal  amplitude  and  time  delay 
measurements . 

The  pulse  is  also  displayed  on  an  oscilloscope.  The  tone  burst 
generator  has  a  synchronized  output,  which  is  connected  to  the  external 
trigger  input  of  the  oscilloscope.  This  causes  the  oscilloscope  trace 
to  trigger  at  the  beginning  of  the  four-cycle  pulse. 

The  acoustic  pulse  generated  in  the  sediment  is  detected  by  the 
receiving  transducer.  The  transducer  converts  the  acoustic  signal  into 
a  voltage  which  is  sent  through  a  bandpass  filter,  and  then  displayed  on 
the  oscilloscope.  The  low  and  high  cutoff  frequencies  of  the  bandpass 
filter  are  100  kHz  and  110  kHz,  respectively.  The  filter  is  used  to 
remove  the  extraneous  frequencies,  introduced  because  the  transmitted 
pulse  is  the  product  of  a  pure  tone  and  the  difference  of  two  step 
functions,  from  the  received  signal.  The  measurements  of  attenuation 
and  phase  velocity  will  then  correspond  to  the  single  frequency  of 
105.0  kHz. 

I  C .  T ransducers 

1 .  Construction 

The  two  transducers  are  nearly  identical  and  can  be  used  as 
either  a  transmitter  or  receiver.  The  transducers  were  built  from  a 
piezoelectric  element  mounted  in  a  stainless  steel  plate.  A  rod  was 
welded  onto  the  plate  so  that  the  plate  and  transducer  element  could  be 
shoved  into  the  sed iment-water  mixture.  A  sketch  of  a  transducer  mount 
is  shown  in  Fig.  2. 

The  two  elements  are  made  of  Channel ite  5500,  manufactured  by 
Channel  Industries,  Inc.  Their  thickness  and  diameter  are  3/16  in.  and 
3/4  in.,  respectively.  Each  element  has  two  resonances.  In  the 

thickness  mode,  transducers  No.  1  and  2  resonate  at  446.4  kHz  and 
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445.2  kHz,  respectively.  In  the  radial  mode,  transducers  No.  1  and  2 
resonate  at  105.5  kHz  and  105.1  kHz,  respecti vely.  In  this  experiment, 
the  transducers  are  used  in  the  radial  mode,  at  105.0  kHz. 

To  build  the  transducers,  a  7/8  in.  hole  was  drilled  through 
the  steel  plate.  The  element  was  mounted  in  the  hole  with  a  layer  of 
corprene  placed  between  it  and  the  steel  plate;  a  layer  of  corprene  was 
also  attached  on  the  backside  of  the  transducer.  A  wire  was  soldered  on 
each  face  of  the  element,  and  the  two  wires  were  placed  in  a  shielded, 
waterproofed  sheath. 

Corprene,  which  has  a  much  lower  acoustical  impedance  than  the 
piezoelectric  element  or  the  sediment-water  mixture,  is  used  on  the 
backside  of  the  element  as  a  pressure  release  material.  Thus  most  of 
the  energy  is  radiated  from  the  front  side  of  the  transducer,  toward  the 
receiver,  thereby  producing  a  higher  signal -to-noise  ratio. 

The  plate  and  the  element  were  designed  to  act  as  a  radiating 
piston  in  a  rigid  baffle.  In  this  way,  the  geometrical  spreading  and 
the  directivity  of  the  transducer  could  be  calculated  and  verified  by 
measurement. 

The  circular  radiating  surface  was  large  enough  to  allow  a 
substantial  transfer  of  energy  from  the  piezoelectric  element  to  the 
sediment.  Since  this  surface  area  was  large  compared  to  the  individual 
sediment  grains,  problems  associated  with  the  coupling  of  the  transducer 
to  the  sediment  were  reduced. 

2 .  Geometrical  Spreading 

Since  the  geometrical  spreading  in  the  nearfield  is  very 
complicated  and  varies  with  distance  from  the  source,  the  transducers 
were  designed  for  use  in  the  farfield.  The  Rayleigh  distance  is  roughly 
the  distance  away  from  the  source  at  which  the  farfield  begins.  At 
distances  greater  than  the  Rayleigh  distance,  the  amplitude  of  the 
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signal  decays  as  1/r,  due  to  the  geometrical  spreading  of  the  signal. 
Here,  r  is  the  distance  from  the  source. 

The  Rayleigh  distance  for  a  circular  piston  in  an  infinite, 
18 

rigid  baffle  is  given  by 


where  a  is  the  radius  of  the  piston,  k=w/c  is  the  wave  number,  and  w  is 
the  angular  frequency.  Note  that  R  is  inversely  proportional  to  the 
phase  velocity  c.  For  these  transducers,  operating  at  a  frequency  of 
105.0  kHz,  this  formula  gives  a  Rayleigh  distance  of  2.0  cm  in  water. 
Spherical  spreading  is  expected  at  separations  of  about  4.0  cm  or 
greater. 

For  a  lossless  medium  in  which  there  is  1/r  spreading,  the 
signal  amplitude  satisfies  the  following  relation. 

E/Er  =  KA/r 

where  E  is  the  signal  amplitude,  Er  is  the  reference  amplitude,  A  is  the 
wavelength,  and  K  is  a  calibration  constant.  Taking  the  logarithm  of 
both  sides  gives 

E  r 

log  r  =  -log  r-  +  log  K 
Lr  A 

Thus,  if  the  signal  amplitude  is  measured  as  a  function  of  separation 

for  separations  greater  than  the  Rayleigh  distance,  the  points  will  fall 
on  a  straight  line  with  a  slope  of  -1  when  plotted  on  log-log  graph 

paper. 

To  experimentally  verify  that  there  is  spherical  spreading  in 

water  at  separations  greater  than  4.0  cm,  the  amplitude  of  the  acoustic 

signal  was  measured  as  the  range  was  varied  from  4.0  cm  to  20.0  cm,  in 


2.0  cm  increments.  A  least  squares  fit  through  the  points  on  a  log-log 
plot  of  range  versus  amplitude  reveals  a  slope  of  m=-0.9914  of  the 
equation 


r  p 

log  £  =  m  log  j-  +  log  B 

If  the  point  corresponding  to  r  =  4.0  cm  is  omitted,  the  slope  is 
m=-1.0061  as  shown  in  Fig.  3.  At  these  distances,  it  can  be  assumed 
that  there  is  no  absorption  in  water  and  the  total  losses  can  be 
attributed  to  geometrical  spreading.  Therefore,  this  test  shows  that 
there  is  spherical  spreading  at  separations  greater  than  4.0  cm. 

Since  the  phase  velocity  is  greater  in  the  sediment-water 
mixture  than  in  water,  the  Rayleigh  distance  is  even  smaller  in  the 
mixture.  Thus  it  is  known  that  the  geometrical  spreading  in  the 
sediment  is  spherical  at  separations  greater  than  4  cm. 

3.  Directivity 


Since  the  transducers  are  directional,  great  care  has  to  be 
taken  in  moving  the  receiving  transducer  to  keep  the  transducers  on  each 
other's  axes  where  the  signal  amplitude  is  maximum.  The  axis  is  the 
imaginary  line  running  through  the  center  of  the  circular  element  and 
perpendicular  to  its  face.  The  alignment  problem  is  lessened  if  the 
radiation  pattern  of  the  transducer  is  sufficiently  wide  that  slight 
errors  in  realigning  the  receiver  exactly  on  the  axis  of  the  transmitter 
will  cause  only  a  slight  error  in  the  on-axis  signal  amplitude. 

By  treating  the  transducer  as  a  circular  piston  in  an 
infinite,  rigid  baffle,  the  half-power  angle  can  be  calculated.  If 
p(0,r,t)  is  the  on-axis  pressure,  then  the  half-power  angle  is  the  angle 
at  which  the  pressure  is  given  by^ 

p(0uP,r,t)  =  A-  P(0,r,t) 

/  2 


For  a  circular  piston,  the  half-power  angle  is 
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s  i  n0|_|  p  =  1.616/(ka) 


In  fresh  water,  at  a  frequency  of  105.0  kHz,  the  calculated  half-power 
angle  for  these  transducers  is  22.6°. 

Figures  4  and  5  show  the  measured  directivity  of  the  two 
transducers  in  fresh  water.  Again,  it  is  seen  that  the  model  of  the 
transducer  as  a  circular  piston  on  a  rigid  baffle  works  well.  The 
directivity  patterns  for  both  transducers  show  that  the  half-power 
angle,  the  angle  at  which  the  pressure  is  3.0  dB  below  the  on-axis 
pressure,  is  about  20°. 

Since  transducer  No.  1  is  used  as  the  receiver  in  this 
experiment ,  and  is  the  one  being  moved  after  each  reading,  it  is  desired 
that  the  response  of  this  transducer  especially  be  constant  for  small 
angles  off  axis.  This  is  so  that  failure  to  exactly  realign  the 
transducers  will  still  give  a  reliable  measurement  of  the  on-axis  value 
of  the  signal  amplitude.  Figure  4  shows  that  the  sound  pressure  level 
for  transducer  No.  1  varies  less  than  0.5  dB  for  7.0°  on  either  side  of 
the  axis.  A  0.5  dB  difference  in  sound  pressure  level  implies  that  the 
actual  pressure  varies  less  than  6.0%  from  the  on-axis  value.  Actually, 
the  SPl  for  transducer  No.  1  is  essentially  constant  for  angles  up  to 
3.0°  off  axis. 

For  mediums  such  as  saturated  sediments,  in  which  the  sound 
speed  is  higher  than  that  of  water,  the  half-power  angle  would  be  even 
greater.  A  typical  compressional  wave  speed  for  a  sediment-water 
mixture  is  1800  m/s. ^  At  the  same  frequency  that  is  being  considered 
(105.0  kHz),  this  means  a  half-power  angle  of  28°.  Thus,  for  small 
angles  off  axis,  the  signal  amplitude  would  vary  even  less  from  the  on- 
axis  value  in  the  mixture  than  it  would  in  water.  In  the  measurements 
taken  in  this  experiment,  care  was  taken  to  position  the  receiver  in  the 
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sediment  with  no  more  than  4  variation  in  the  angle  off  axis.  Thus,  we 
have  confidence  that  the  on-axis  signal  amplitudes  were  measured 
accurately. 

D.  Sediment  Preparation 

The  sediment  is  made  up  of  MS-MH  size  glass  beads  manufactured  by 
Ferro  Corporation,  Jackson,  Mississippi.  The  grain  size  distribution  is 
shown  in  Fig.  6. ^ 

The  sediment-water  mixture  was  prepared  in  a  5  gal  plastic  bucket 
in  the  following  manner.  The  bucket  is  approximately  35  cm  high  and 

29  cm  in  diameter.  The  liquid  in  the  sediment-water  mixture  is 
circulated  by  a  Masterflex  tubing  pump.  Before  the  mixture  was  put  in 
the  bucket,  the  outflow  tube  was  attached  to  the  bottom  of  the  bucket 
with  duct  tape.  The  bucket  was  then  partially  filled  with  de-ionized 
water  and  the  input  hose  was  put  in  the  water.  The  pump  was  turned  on 
so  that  all  of  the  air  was  driven  from  the  tube.  A  capful  of  BIO¬ 
CLEAN  III  algaecide  was  added  to  the  water  to  retard  the  growth  of 

micro-organisms . 

The  glass  shot  was  then  poured  slowly  into  the  water  until  the 

sediment-water  mixture  was  24  cm  deep.  Pouring  was  halted  at  intervals 
and  the  mixture  kneaded  like  bread  dough  to  remove  air  trapped  in  the 

beads.  A  layer  of  water,  approximately  6  cm  thick,  was  left  on  top  of 

the  sediment  mixture  and  the  pump  input  hose  was  suspended  in  the  layer, 
as  shown  in  Fig.  7. 

To  remove  any  air  remaining  in  the  mixture,  the  mixture  was 

vibrated  in  excess  of  6  hours.  The  mixture  was  left  undisturbed 

overnight.  The  next  day,  the  transducers  were  shoved  into  the  mixture 
so  that  the  centers  of  the  elements  were  12  cm  into  the  sediment. 


To  increase  the  glycerine  concentration  of  the  liquid,  glycerine 
was  poured  into  the  liquid  on  top  of  the  mixture  and  the  pump  was  turned 
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FIGURE  6 

GRAIN  SIZE  DISTRIBUTION  FOR  MS-MH  SIZE  BEADS 


FIGURE  7 

CIRCULATING  THE  PORE  FLUID 


on  and  run  until  the  liquid  in  the  bucket  was  well  mixed.  A  capful  of 

algaecide  was  added  while  the  liquid  was  being  circulated  to  further 

retard  the  growth  of  micro-organisms.  The  specific  gravity  of  the 

liquid  was  measured  several  times  while  the  liquid  was  being  circulated. 

The  liquid  was  considered  well  mixed  after  the  values  of  specific 

gravity  had  stabilized.  This  usually  occurred  after  about  3  hours.  The 

specific  gravity  was  measured  using  a  hydrometer  and  the  value  of 

viscosity  was  read  from  a  table  containing  the  values  of  viscosity  for 
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varying  concentrations  of  glycerine-water  solutions. 

E .  Measurements  and  Preliminary  Results 
1 .  Phase  Velocity 

The  oscilloscope  trace  was  triggered  at  the  beginning  of  the 
transmitted  four-cycle  pulse.  The  arrival  time  to  the  fourth  positive 
peak  on  the  received  signal  was  measured  with  a  digital  readout  on  the 

oscilloscope.  This  was  first  done  at  a  transducer  separation  of 

approximately  20  cm  and  then  at  2.0  cm  increments  down  to  a  separation 
of  6  cm. 

If  these  arrival  times  were  then  plotted  on  a  time  versus  distance 
graph,  with  time  as  the  horizontal  axis,  the  points  should  fall  on  a 
straight  line  parallel  to  the  characteristic  of  the  wavelet.  The  slope 
of  this  line  is  the  phase  velocity.  Therefore,  to  find  the  phase 
velocity  a  least  squares  routine  is  used  to  fit  the  best  straight  line 
through  these  points.  Since  the  line  is  parallel  to  a  characteristic, 
the  equation  of  the  line  is  of  the  form 

r  =  ct  +  B  , 


where  c  is  the  phase  velocity  and  B  is  the  intercept  of  the  r  axis.  If 
the  time  were  measured  from  the  fourth  positive  peak  on  the  transmitted 
pulse  to  the  same  peak  on  the  received  pulse,  B  would  ideally  be  zero 
(for  a  nondispersi ve  medium)  and  the  least  squares  line  would  be  the 


same  as  the  characteristic.  But  since  time  is  measured  from  the 
beginning  of  the  triggered  pulse,  and  since  there  are  additional  time 
delays  associated  with  the  circuitry  and  transducers,  B  turns  out  to  be 
negative. 

This  method  was  tested  using  two  different  size  beads  with  fresh 
water  as  the  pore  fluid:  the  MS-MH  size,  which  was  used  in  the 

experiment,  and  the  MS-L  size  beads,  which  are  a  smaller  size.  The 
results  of  these  preliminary  tests  are  shown  in  Figs.  8  and  9.  The 
linear  correlation  coefficient  is  a  number  between  0  and  1  that  gives  an 
indication  of  how  well  the  data  points  fit  a  straight  line.  A 

coefficient  of  1  represents  perfect  linear  correlation  and  a  coefficient 
of  0  represents  no  correlation.  It  can  be  seen  in  both  of  these  figures 
that  the  correlation  between  the  data  and  the  straight  line  is  very 

good.  The  good  correlation  indicates  that  the  disturbances  in  the 
sediment,  caused  by  moving  the  transducer  to  each  separation,  has  a 
negligible  effect  on  the  measurements. 

In  order  to  check  the  validity  of  our  procedure,  we  measured  the 
phase  velocity  in  water  at  room  temperature,  approximately  25°C,  six 
times.  Each  of  these  results  was  between  1492  m/s  and  1508  m/s.  The 
average  of  these  measurements  was  1501  m/s.  By  linear  interpolation  of 
tabulated  values  of  sound  speed  as  a  function  of  temperature,  we  found 
that  the  speed  of  sound  at  25°C  in  fresh  water  is  1496  m/s.  The  average 
then  has  an  error  of  only  0.33%  while  the  extreme  values  have  errors  of 
-0.3%  and  0.8%. 

In  another  check,  we  compared  our  results  to  results  obtained  by 
Bell ,  ^  Shirley,^  and  Hovem  and  Ingram.^  We  measured  the  phase 
velocity  three  separate  times  in  MS-MH  beads  with  fresh  water  as  the 
pore  fluid.  The  three  results  are  shown  in  Table  I,  along  with  the 
other  results  mentioned.  Comparison  with  these  other  results  indicates 
that  our  results  are  valid.  But  the  comparison  also  suggests  that  the 

velocity  is  strongly  dependent  on  the  depth  of  the  transducers  in  the 

sediment. 
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BEAD  TYPE:  MS-MH 

FREQUENCY:  105  kHz  VISCOSITY:  0.893  cP 


TIME  -  /is 

65.88 

77.63 

88.98 

100.12 

112.07 

121.38 

133.23 

142.31 


SEPARATION  - 

5.65 

7.65 

9.65 

11.65 

13.65 

15.65 

17.65 

19.65 


VELOCITY:  1820  m/s 

STANDARD  DEVIATION  OF  VELOCITY:  24  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9995 
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FIGURE  8 

PHASE  VELOCITY  IN  MS-MH  BEADS 
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BEAD  TYPE:  MS-L 

FREQUENCY:  105.0  kHz  VISCOSITY:  0.893  cP 


l/IE  —  us 

SEPARATION  -  cm 

63.62 

5.65 

73.39 

7.65 

83.49 

9.65 

95.51 

11.65 

107.15 

13.65 

118.60 

15.65 

131.73 

17.65 

138.98 

19.65 

VELOCITY:  1791  m/s 

STANDARD  DEVIATION  OF  VELOCITY:  34  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9989 


FIGURE  9 

PHASE  VELOCITY  IN  MS-L  BEADS 
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TABLE  I 

COMPARISON  OF  VELOCITY  RESULTS  IN  MS-MH  AND  MS-L  SIZE  BEADS 


Frequency 

(kHz) 


Bell 

120.0 

Shirley 

114.0 

Hovem  &  Ingram 

100.0 

Present  Study 

105.0 

Depth  Velocity 

(cm)  (m/s) 


MS-MH 

MS 

9 

1809 

13 

10 

1808 

20 

1922 

12 

1820 

17 

The  comparison  in  the  MS-L  size  beads  is  not  nearly  so  good.  Our 
first  trial  measurement  in  sediments  was  made  in  the  MS-L  size  beads. 
Since  the  standard  deviation  accompanying  this  measurement  is  so  large 
(34  m/s)  and  since  in  all  subsequent  measurements  the  standard  deviation 
was  substantially  reduced,  the  discrepancy  is  attributed  to  experimental 
error  and  carelessness  in  precisely  determining  the  separation.  As  the 
experiment  progressed,  the  technique  improved  and  the  error  decreased. 

Obviously,  sound  speed  measurements  in  sediments  are  not  nearly  as 
repeatable  as  measurements  made  in  water.  One  reason  for  this  is  the 
uncertainty  in  the  separation  measurement.  In  the  sediment,  this 
uncertainty  is  estimated  to  be  0.1  cm.  At  a  separation  of  20  cm,  this 
represents  a  relative  uncertainty  of  0.5%,  while  at  6  cm  the  uncertainty 
is  1.7%.  The  uncertainty  in  the  time  delay  measurement  is  negligible 
compared  to  this.  Relative  uncertainties  of  0.5%  and  1.7%  in  the 
separation  represent  absolute  uncertainties  in  the  velocity  of  10  m/s 
and  30  m/s,  respecti vely.  The  standard  deviations  shown  in  Figs.  8 
and  9  were  derived  from  the  least  squares  method.  The  assumptions  made 
in  these  calculations  were  that  the  absolute  uncertainties  in  the 

vertical  axis,  or  separations,  are  equal,  and  the  uncertainties  in  the 
horizontal  axis  are  negligible.  Both  of  these  assumptions  are 
appropriate,  and  the  standard  deviations  calculated  from  this  method  are 
consistent  with  uncertainties  in  separation  and  hence  are  a  realistic 
measure  of  the  uncertainty  in  the  velocity. 

The  standard  deviations  of  the  velocities  were  less  than  15  m/s 

after  the  first  trial  measurements.  This  indicates  that  as  the 
experiment  progressed,  the  technique  improved  and  the  uncertainty 

decreased. 

2 .  Attenuation 

In  the  sediment,  the  losses  are  due  to  both  geometrical 
spreading  and  absorption  of  the  medium.  Therefore,  to  calculate  the 

attenuation  in  the  medium,  the  geometrical  spreading  must  be  taken  into 


account.  As  discussed  in  Section  II. C. 2,  the  transducers  were  designed 
so  that  the  Rayleigh  distance  would  be  relatively  small  compared  to  the 
sediment  sample.  At  distances  greater  than  the  Rayleigh  distance,  the 
signal  amplitude  decays  as  1/r  due  to  the  spherical  spreading,  where  r 
is  the  separation  between  the  transducers.  This  distance  was  calculated 
for  both  water  and  sediment.  The  calculation  for  water  was  then 
verified  by  acoustic  measurements.  These  measurements  show  that  the 
geometrical  spreading  in  the  sediment  was  indeed  spherical  at  the 
distances  at  which  the  attenuation  measurements  were  taken.  Using  the 
spherical  spreading,  the  losses  due  to  absorption  in  the  sediment  can  be 
readily  determined.  A  detailed  discussion  of  this  procedure  follows. 

The  measurements  establishing  the  1/r  decay  of  the  signal 
amplitude  due  to  geometrical  spreading  has  been  discussed  in 
Section  II. C. 2.  The  absorption  due  to  the  medium  can  be  determined  by 
measuring  the  amplitude  of  the  signal  and  assuming  that  the  amplitude 
varies  as 


r  rv 

E  =  ?  e 


Here,  K  is  a  calibration  constant  with  units  of  volts  x  length  and  a  is 
the  attenuation  in  nepers  per  unit  length.  The  1/r  factor  in  Eq.  (1) 
takes  into  account  the  spherical  spreading  and  the  exponential  term 
takes  into  account  the  absorption  of  the  medium.  It  is  also  assumed 
that  the  coupling  between  the  sediment  and  the  transducer  is  the  same  at 
each  separation.  Thus,  the  decay  in  the  amplitude  of  the  signal  is 
attributed  only  to  the  increase  in  separation. 

The  attenuation  is  determined  by  comparing  the  amplitudes 
measured  at  various  separations.  To  do  this,  one  of  the  amplitude 
measurements  and  the  separation  at  which  it  is  measured  are  designated 
as  reference  values  Er  and  r  ,  respectively.  In  this  experiment,  rr  is 
the  smallest  separation  at  which  measurements  are  taken  and  E  =E  at 


Therefore,  the  following  relation  holds. 


f-  in 

'•  =  ve 


From  this  and  Eq.  (1)  it  follows  that: 

a  (r-rf)  =  in  •  (2) 

Therefore,  if  Jin  ^(Errr)/Erj  versus  (r-rr)  is  plotted,  the  points 
should  fall  on  a  straight  line  and  the  slope  of  the  line  would  be  the 

attenuation.  To  find  the  slope  of  this  line,  and  thereby  the 

attenuation,  the  method  of  least  squares  is  used. 

This  technique  was  tested  using  two  different  size  beads  with 
fresh  water  as  the  pore  fluid:  the  MS-MH  size,  which  was  used  in  the 

experiment,  and  the  MS-L  size,  which  are  smaller  beads.  The  results  of 

these  preliminary  tests  are  shown  in  Figs.  10  and  11.  As  in  the  phase 
velocity  case,  a  linear  correlation  coefficient  of  near  1  indicates  a 
very  close  correlation  between  the  data  and  the  straight  line. 

The  standard  deviation  can  be  calculated  for  this  method, 
assuming  that  the  uncertainties  in  all  of  the  data  points  are  equal  and 
that  there  is  no  error  in  the  values  along  the  horizontal  axis.  Because 
these  assumptions  are  false,  the  standard  deviation  calculated  from  the 
least  squares  method  is  too  low.  But  the  attenuation  can  be  calculated 
by  other  means,  and  from  these  calculations  a  realistic  estimate  of  the 
error  can  be  made. 

A  physical  interpretation  of  the  sources  of  error  can  be 
acquired  by  calculating  the  attenuation  directly,  using  Eq.  (2)  with  the 
data  at  any  two  separations.  If  the  data  for  MS-MH  beads  corresponding 
to  separations  of  5.65  cm  and  19.65  cm  are  used,  as  shown  in  Fig.  11, 
the  attenuation  turns  out  to  be  2.28  Np/m.  We  can  calculate  a 
Pythagorean  estimate  of  the  uncertainty  in  the  attenuation  by  using  the 
relation2^ 


BEAD  TYPE:  MS-L 


FREQUENCY:  105.0  kHz 

VISCOSITY: 

AMPLITUDE  -V 

SEPARATION 

0.2800 

5.65 

0.1970 

7.65 

0.1430 

9.65 

0.1060 

11.65 

0.0820 

13.65 

0.0652 

15.65 

0.0516 

17.65 

0.0430 

19.65 

ATTENUATION:  4.71  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.14  Np/m 
LINEAR  CORRELATION  COEFFICIENT:  0.9974 

E  =  AMPLITUDE  r  =  SEPARATION 


FIGURE  10 

ATTENUATION  IN  MS-L  BEADS 
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BEAD  TYPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  0.893  cP 


AMPLITUDE  -  V 

0.3720 

0.2690 

0.2004 

0.1570 

0.1300 

0.1074 

0.0892 

0.0777 


SEPARATION  -  cm 

5.65 

7.65 

9.65 

11.65 

13.65 

15.65 

17.65 

19.65 


ATTENUATION:  2.40  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.09  Np/m 
LINEAR  CORRELATION  COEFFICIENT:  0.9962 


FIGURE  11 

ATTENUATION  IN  MS-MH  BEADS 
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a- (r-rr)2  =  a2(a2  +  a2  a/r  +  l/r2)ref  +  a2(a2  +  2a /r  +  1/r2)  ,  (3) 

+  {aE/E)?ef  +  (0E/E)2 

where  ac  and  o  denote  the  uncertainties,  or  standard  deviations,  in  the 
c  r 

amplitude  and  separation,  respecti vely ,  and  where  the  quantities  marked 
"ref"  signify  that  the  separations  and  amplitudes  correspond  to  the 
separation  of  5.65  cm.  To  perform  this  calculation,  the  uncertainties 
in  the  separation  and  the  amplitude,  and  cr^,  must  first  be 

determined.  The  uncertainty  in  the  separation  is  estimated  to  be 
0.1  cm. 


The  uncertainty  in  the  amplitude  has  two  sources.  The  first 

one  is  instrumental;  that  is,  the  oscilloscope  can  be  read  accurately  to 

1/4  of  a  gradation.  Because  the  scale  is  changed  at  different 
separations,  since  the  amplitude  grows  with  decreasing  r,  this  type  of 
uncertainty  increases  with  decreasing  separation. 

The  second  source  of  uncertainty  in  E  is  due  to  the 
uncertainty  in  r.  Since  the  amplitude  is  a  function  of  r,  as  was 
assumed  in  Eq.  (1),  any  uncertainty  in  r  produces  an  uncertainty  in  E. 

This  type  of  error  can  be  estimated  from  the  relation 

a  =  a  ~  =  0  E(ar  +  I) 

E  r  a  r  r  v  r' 

By  calculating  this  uncertainty  and  comparing  the  results  with  the 

instrumental  uncertainty  at  each  separation,  it  can  be  seen  that  over  W; 

80%  of  the  uncertainty  in  E  is  due  to  the  uncertainty  in  r.  Using 
Eq.  (3)  and  the  values  in  the  rows  in  Table  II  corresponding  to 
r  =  19.65  cm  and  r  =  5.65  cm,  the  uncertainty  in  the  direct  calculation 
is  found  to  be  0.522  Np/m.  Note  from  Eq.  (3)  that  the  uncertainty  a.  is 
inversely  proportional  to  the  difference  between  the  two  separations 
being  compared.  Thus,  one  way  to  improve  the  measurements  would  be  to 
increase  the  distances  over  which  the  measurements  are  taken. 
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Using  each  pair  of  data  shown  in  Fig.  11,  the  attenuation  can 

be  calculated  seven  times  in  this  way,  and  the  values  designated  as  . 

These  seven  values  of  attenuation  can  then  be  weighted  with  their 

20 

corresponding  uncertainties  a.,  and  averaged  as  shown  in  Table  II.  The 
resulting  average  is  a  better  estimate  of  the  attenuation  than  any  of 
the  individual  attenuations,  and  the  standard  deviation  calculated  from 
this  method  is  a  better  estimate  of  the  uncertainty.  Although  this 
value  of  a  differs  by  5%  from  the  least  squares  value,  this  standard 
deviation  is  a  more  realistic  estimate  of  the  uncertainty. 

In  calculating  the  attenuation  using  this  method,  the  error  in 
the  reference  separation  and  amplitude  is  repeated  in  each  of  the  seven 
calculations.  Thus,  if  that  measurement  is  particularly  bad,  the  error 
in  it  will  be  compounded  in  the  final  result.  To  avoid  this  type  of 
error,  the  attenuation  can  be  calculated  with  any  two  pairs  of  data,  E 
and  r.  The  values  corresponding  to  r  =  5.65  cm  do  not  necessarily  need 
to  be  the  reference  separation  and  amplitude.  Since  there  are  eight 
pairs  of  data,  there  are  28  unique  ways  to  calculate  using  Eq.  (2). 
The  unweighted  mean  of  these  28  attenuations  and  the  standard  deviation 
of  this  mean  can  easily  be  calculated.  For  the  MS-MH  beads  of  Fig.  11, 
this  mean  is  a  =  2.38  Np/m  and  the  standard  deviation  isa=  0.47  Np/m. 
This  value  of  the  attenuation  differs  by  only  0.6%  from  the  least 
squares  value.  But  this  standard  deviation  is  larger  than  the  one 
calculated  from  the  weighted  averages.  This  may  signify  that  the 
estimate  of  the  uncertainty  in  r  used  in  that  calculation,  =  0.1  cm, 
was  too  low.  In  later  calculations,  however,  there  is  better  agreement 
between  these  two  calculations  of  the  standard  deviation.  This  suggests 
that  the  separation  was  determined  more  precisely  as  the  experiment 
progressed.  Another  reason  for  the  discrepancy  is  that  all  of  the 
attenuations  are  weighted  equally  in  calculating  the  mean.  But  as  can 
be  seen  from  Eq.  (3),  the  uncertainty  decreases  as  the  difference 
between  the  two  separations  increases.  Thus,  the  standard  deviation 
calculated  by  this  method  represents  an  upper  bound  of  the  actual 
uncertainty  in  the  attenuation. 
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The  attenuation  was  calculated  for  each  value  of  viscosity 

using  the  method  of  least  squares  and  the  mean  of  all  combinations.  In 

each  case,  the  two  values  of  attenuation  differed  by  less  than  1%.  But 

we  feel  that  the  standard  deviation  associated  with  the  mean  of  all 
combinations  was  more  realistic. 

It  can  be  noted  that  the  uncertainty  in  the  separation  is  by 
far  the  biggest  source  of  error.  Therefore,  to  improve  these 

measurements,  a  more  accurate  method  for  determining  the  separation  of 
the  transducers  would  have  to  be  devised.  Toward  this  end,  the 
separations  were  corrected  using  the  velocity  results.  The  attenuations 
were  then  calculated  by  the  mean  of  all  combinations  with  these 
corrected  separations.  These  attenuations  were  within  \%  of  the 
attenuations  calculated  by  the  other  methods,  except  for  the  case  of  the 
MS-MH  beads  of  Fig.  11  where  the  difference  was  2.5%.  In  all  but  a  few 
cases,  the  standard  deviation  was  substantially  reduced.  Thus,  the 
attenuation  was  calculated  by  the  mean  of  all  combinations.  In  all  but 
those  few  cases  where  the  standard  deviation  was  not  reduced,  the 
separations  were  corrected  using  the  velocity  results.  Reducing  the 
data  in  this  way,  we  are  confident  that  the  attenuations  were  determined 
within  3  dB/m,  except  for  a  few  cases. 

3.  Comoarison  of  Results:  Attenuation 

To  determine  the  validity  of  the  attenuation  results,  the 
measurements  in  the  water  saturated  beads,  with  no  glycerine  added,  were 
compared  to  the  results  of  David  Bell .  ^  His  attenuation  measurements 
are  given  in  decibels  per  wavelength.  Converting  his  results  to  nepers 

per  meter  at  105.0  kHz  gives  the  values  5.45  Np/m  and  3.13  Np/m  for  MS-L 

3 

and  MS-MH  size  beads,  respectively. 

The  discrepancy  in  attenuation  between  Bell's  values  and  those 
presented  in  Figs.  10  and  11  is  0.74  Np/m  for  MS-L  and  0.73  Np/m  for 

MS-MS.  Since  this  difference  is  too  large  to  be  attributed  to 

experimental  error,  an  effort  was  made  to  explain  it. 


In  order  to  account  for  the  geometrical  spreading  of  the 

transducers,  Bell  measured  the  amplitude  decay  in  water  at  120  kHz  and 
the  amplitude  decay  in  the  sediment  at  the  same  frequency.  He  assumed 

that  the  only  losses  in  the  water  were  due  to  geometrical  spreading  and 

that  the  losses  in  the  sediment  were  due  to  geometrical  spreading  and 

the  attenuation.  He  then  assumed  that  the  losses  due  to  geometrical 
spreading  in  the  sediment  were  the  same  as  in  the  water  at  that  same 
frequency.  Therefore,  to  calculate  the  attenuation  in  the  medium  he 

subtracted  the  spreading  losses,  which  he  determined  from  his 

measurements  in  water,  from  the  total  losses. 

In  doing  this,  he  assumed  that  the  geometrical  spreading  in 

the  water  was  the  same  as  in  the  sediment  for  measurements  taken  at  the 
same  frequency.  But  the  Rayleigh  distance,  which  is  a  measure  of  the 
distance  at  which  spherical  spreading  starts,  is  inversely  proportional 
to  the  sound  speed.  Since  the  sound  speed  is  greater  in  the  sediment 
than  in  the  water,  the  Rayleigh  distance  is  smaller  in  the  sediment  than 
in  the  water.  Also,  the  distances  at  which  he  took  his  measurements 
were  less  than  the  Rayleigh  distance  for  his  transducers,  so  that  in 

water  his  amplitude  did  not  vary  as  1/r.  Therefore,  the  losses  due  to 
geometrical  spreading  were  less  in  the  water  than  in  the  sediment.  The 
extra  geometrical  spreading  that  was  occurring  in  the  sediment  was 
attributed  to  absorption  in  the  sediment  resulting  in  a  value  of 
attenuation  that  was  too  large. 

We  also  compared  our  measurements  in  MS-MH  beads  with  those 
made  by  Shirley  ^  and  Hovem  and  Ingram.**  These  results  are  summarized 
in  Table  III.  Shirley's  results,  like  Bell's,  are  much  higher  than 
ours.  We  believe  that  he  measured  the  attenuation  in  the  same  way  that 
Bell  did.  Thus  his  results  differ  from  ours  for  the  same  reason  that 
Bell's  do. 


4 

Hovem  and  Ingram's  results  at  100  kHz  also  differ 
significantly  from  ours.  If  their  results  at  only  this  one  frequency 
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TABLE  III 

COMPARISON  OF  ATTENUATION  RESULTS  IN  MS-MH  SIZE  BEADS 


Frequency 

Depth 

Attenuation 

(kHz) 

(cm) 

(dB/m) 

Present  Study 

- 5 — 1 — 

105.0 

12.0 

20.85 

20.74 

21.61 

Bell 

120.0 

9.0 

27.10 

Shirley 

114.0 

10.0 

30.91 

Hovem  X  Ingram 

60.0 

20.0 

24.86 

80.0 

20.0 

17.63 

100.0 

20.0 

13.20 

150.0 

20.0 

19.65 
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are  considered,  one  might  be  led  to  the  conclusion  that  the  measured 
attenuation  is  very  much  dependent  on  the  depth  of  the  transducers  in 
the  sediment.  But  if  their  results  at  other  frequencies  in  that  range 
are  considered,  it  can  be  seen  that  their  data  exhibit  a  significant 
amount  of  dispersion.  It  is  difficult  to  say  whether  this  dispersion  is 
due  to  actual  physical  processes  or  to  scatter  in  the  data.  Other 
frequency  dependent  measurements  need  to  be  taken  to  investigate  this. 

There  were  also  several  assumptions  made  in  our  measurements 
which  could  affect  the  validity  of  the  attenuation  results.  One 
assumption  was  that  the  coupling  between  the  transducer  and  the  sediment 
was  the  same  in  each  case.  Since  this  is  not  a  trivial  assumption,  an 
effort  was  made  to  justify  it.  Each  time  the  transducer  was  moved,  the 
sediment-water  mixture  was  vibrated  gently  for  approximately  30  minutes. 
The  mixture  was  then  left  undisturbed  for  an  additional  half  hour  before 
another  measurement  was  taken.  Four  or  five  more  readings  were  then 
taken  over  the  next  3-4  hours.  Since  the  voltage  and  travel  time 
readings  never  changed  over  this  period  of  time,  except  for  small 
variations  attributed  to  experimental  uncertainty,  it  was  assumed  that 
the  coupling  between  the  transducer  and  sediment  did  not  vary. 


Summan 


The  measurement  techniques  described  in  this  report  are  refinements 
of  techniques  developed  and  used  previously.^’ ^  The  main  effort, 
thus  far,  has  been  to  design  a  valid  technique  for  measuring  the  phase 
velocity  and  attenuation  of  compressional  waves  in  sediments.  The  main 
improvement  over  the  earlier  work  has  resulted  from  redesigning  the 
compressional  wave  transducers  so  that  the  geometrical  spreading  and 
directivity  can  be  predicted.  With  this  improvement,  spherical 
spreading  is  established  at  the  separations  at  which  these  measurements 
are  taken.  Thus  the  losses  due  to  geometrical  spreading  are  accurately 
accounted  for  when  analyzing  the  data. 


III.  DISCUSSION 


The  purpose  of  the  work  presented  here  has  been  to  design  an 
experiment  that  wi 1 1  examine  the  dependence  of  compressional  wave 
propagation  in  fluid  saturated  sands  on  the  viscosity  of  the  pore  fluid. 
Since  this  dependence  is  predicted  by  the  Biot-Stoll  equations,  these 
measurements  should  be  a  critical  test  of  the  theory. 

1C 

Experiments  of  this  type  have  already  been  performed  by  Shirley, 
and  later  by  ElliotJ®  However,  due  to  the  paramount  difficulties 

encountered  when  taking  measurements  in  sediments,  there  was  too  much 
scatter  in  their  data  to  be  able  to  make  any  definite  conclusion  about 
their  results. 

After  making  a  few  modifications  in  their  technique,  new 
measurements  have  been  taken  and  have  been  shown  to  be  repeatable.  The 
main  modification  has  been  the  redesign  of  the  compressional  wave 
transducer.  The  transducers  were  designed  so  that  the  Rayleigh  distance 
would  be  relatively  short  compared  to  the  size  of  the  sediment  sample. 
Thus,  spherical  spreading  was  established  at  the  distances  at  which  the 
measurements  were  taken.  The  spreading  losses  could  then  be  taken  into 
account  when  analyzing  the  data  so  that  the  attenuation  could  be 
determined. 

An  important  by-product  of  the  new  design  was  the  wide  beam  pattern 
of  the  new  transducer,  which  was  fairly  wide.  As  a  result,  the  sound 
pressure  level  is  essentially  constant  for  small  angles  off  the  axis  of 
the  transducer.  This  enables  the  transducer  to  be  repositioned  at  a  new 
separation  without  the  amplitude  being  affected  by  the  directivity  of 
the  transducer.  Thus,  we  were  confident  that  the  decay  in  amplitude  was 
caused  only  by  increasing  the  separation  of  the  transducers. 

At  this  point,  having  confidence  in  the  measurement  techniques,  the 
phase  velocity  and  attenuation  as  functions  of  viscosity  can  be 
determined.  These  measurements,  once  taken,  can  then  be  compared  to  the 
predictions  of  the  Biot-Stoll  theory. 
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The  University  of  Hawaii 
2525  Correa  Road 
Honolulu,  HI  96822 
Attn:  L.  N.  Frazer 

Director 

North  Atlantic  Treaty  Organization 
SACLANT  ASW  Research  Centre 
APO  New  York  09019 
Attn:  T.  Akal 

F.  Jensen 

Defence  Research  Establishment  Pacific 
FMO  Victoria,  BC 
VOS  130  CANADA 
Attn:  R.  Chapman 
J.  M.  Ozard 
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Defence  Research  Establishment  Atlantic 

9  Grove  Street 
P.  0.  Box  1012 
Dartmouth,  NS 
CANADA 

Attn:  D.  Chapman 

The  Catholic  University  of  America 
6220  Michigan  Avenue,  NE 
Washington,  DC  20017 
Attn:  H.  M.  Uberall 

Lamont-Doherty  Geological  Observatory 
Columbia  University 
Palisades,  NY  10964 
Attn:  G.  Bryan 

R.  D.  Stoll 

Department  of  Civil  and  Ocean  Engineering 
The  University  of  Rhode  Island 
Kingston,  RI  02881 
Attn:  A.  J.  Silva 

University  College  of  North  Wales 
Marine  Science  Laboratories 
Menai  Bridge 
Anglesey,  LL59  5EY,  U.K. 

Attn:  D.  Taylor  Smith 

Institute  of  Oceanographic  Sciences 
Godaiming 

Surrey  GU8  5UB,  U.K. 

Attn:  P.  Schultheiss 

Department  of  Civil  Engineering 
The  University  of  Texas  at  Austin 
Austin,  TX  78712 
Attn.  K.  Stokoe 

Southwest  Research  Institute 
P.  0.  Drawer  28510 
San  Antonio,  TX  78284 
Attn:  D.  J.  Shirley 

Rosenteil  School  of  Marine  and 
Atomospberic  Science 
The  University  of  Miami 

10  Rickenbacker  Causeway 
Miami,  FL  33149 

Attn:  H.  DeFarrari 
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Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Johns  Hopkins  Road 
Laurel,  MD  20810 
Attn:  J.  Lombardo 
R.  Henrick 
A.  Boyles 

Department  of  Ocean  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 
Attn:  I.  Dyer 

G.  Duckworth 
A.  Baggeroer 

The  University  of  Miami 
10  Rickenbacker  Causeway 
Miami,  FL  33149 
Attn:  F.  Tappert 

Department  of  Physics 
The  University  of  Rhode  Island 
Kingston,  RI  02881 
Attn:  C.  Kaufman 

Department  of  Electrical  Engineering 
Polytechnic  Institute  of  New  York 
Farmingdale,  NY  11735 
Attn:  L.  B.  Felsen 

I.  Tolstoy 

Knockvennie,  Castle  Douglas 
S.  W.  SCOTLAND 
GREAT  BRITAIN 

Department  of  Geology 
The  University  of  Texas  at  Austin 
Austin,  TX  78712 
Attn:  C.  Wilson 
M.  Backus 

Department  of  Electrical  Engineering 
The  University  of  Texas  at  Austin 
Austin,  TX  78712 
Attn:  E.  Hixson 

Department  of  Physics 
The  University  of  Auckland 
Private  Bag,  Auckland 
NEW  ZEALAND 

Attn:  A.  C.  Kibblewhite 
C.  T.  Tindle 
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Chinhae  Research  Laboratory 

P.  0.  Box  18 

Chinhae,  Kyeong  Nam 

KOREA 

Attn:  Jungyul  Na 

91 

Department  of  Oceanography 

Texas  A&M  University 

College  Station,  TX  77843 

Attn:  A.  Anderson 

92 

Department  of  Geological  Oceanography 
Texas  A&M  University 

College  Station,  TX  77843 

Attn:  W.  Bryant 

93 

i-95 

Sandia  National  Laboratories 

Division  7552 

Albuquerque,  NM  87185 

Attn:  W.  Shurteleff 

D.  Costley 

96 

A.  Marc  Bedford,  ARL:UT 

97 

H.  Boehme,  ARL:UT 

98 

Robert  A.  Koch,  ARL:UT 

99 

Stephen  K.  Mitchell,  ARL:UT 

100 

T.  G.  Muir,  ARL:UT 

101 

Clark  S.  Penrod,  ARL:UT 

102 

Morris  Stern,  ARL:UT 

103 

Paul  J.  Vidmar,  ARL:UT 

104 

Gary  R.  Wilson,  ARL:UT 

105 

Ching  H.  Yew,  ARL:UT 
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Library,  ARL:UT 

47 


